Evaluating the Performance of the SCOPE Model in Simulating Canopy Solar-Induced Chlorophyll Fluorescence

The SCOPE (soil canopy observation of photochemistry and energy fluxes) model has been widely used to interpret solar-induced chlorophyll fluorescence (SIF) and investigate the SIF-photosynthesis links at different temporal and spatial scales in recent years. In the SCOPE model, the fluorescence quantum efficiency in dark-adapted conditions (FQE) for Photosystem II (fqe2) and Photosystem I (fqe1) were two key parameters of SIF emission, which have always been parameterized as fixed values derived from laboratory measurements. To date, only a few studies have focused on evaluating the SCOPE model for SIF interpretation, and the variation of FQE values in the field remains controversial. In this study, the accuracy of the SCOPE model to simulate the canopy SIF was investigated using diurnal experiments on winter wheat. First, ten diurnal experiments were conducted on winter wheat, and the canopy SIF emissions and the SCOPE model’s input parameters were directly measured or indirectly retrieved from the spectral radiances, gross primary productivity (GPP) data, and meteorological records. Second, the SCOPE-simulated SIF emissions with fixed FQE values were evaluated using the observed canopy SIF data. The results show that the SCOPE model can reliably interpret the diurnal cycles of SIF variation and provide acceptable results of SIF simulations at the O2-B (SIFB) and O2-A (SIFA) bands with RRMSEs of 24.35% and 23.67%, respectively. However, the SCOPE-simulated SIFB and SIFA still contained large systematical deviations at some growth stages of wheat, and the seasonal cycles of the ratio between SIFB and SIFA (SIFA/SIFB) cannot be credibly reproduced. Finally, the SCOPE-simulated SIF emissions with variable FQE values were evaluated using the observed canopy SIF data. The simulating accuracy of SIFB and SIFA can be improved greatly using variable FQE values, and the SCOPE simulations track well with the seasonal SIFA/SIFB values with an RRMSE of 20.63%. The results indicated a clear seasonal pattern of FQE values for unbiased SIF simulation: from the erecting to the flowering stage of wheat, the ratio of fqe1 to fqe2 (fqe1/fqe2) gradually increased from 0.05–0.1 to 0.3–0.5, while the fqe2 value decreased from 0.013 to 0.007. Our quantitative results of the model assessment and the FQE adjustment support the use of the SCOPE model as a powerful tool for interpreting the SIF emissions and can serve as a significant reference for future applications of the SCOPE model

Upscaling Solar-Induced Chlorophyll Fluorescence from an Instantaneous to Daily Scale Gives an Improved Estimation of the Gross Primary Productivity

Solar-induced chlorophyll fluorescence (SIF) is closely linked to the photosynthesis of plants and has the potential to estimate gross primary production (GPP) at different temporal and spatial scales. However, remotely sensed SIF at a ground or space level is usually instantaneous, which cannot represent the daily total SIF. The temporal mismatch between instantaneous SIF (SIFinst) and daily GPP (GPPdaily) impacts their correlation across space and time. Previous studies have upscaled SIFinst to the daily scale based on the diurnal cycle in the cosine of the solar zenith angle ( cos ( SZA ) ) to correct the effects of latitude and length of the day on the variations in the SIF-GPP correlation. However, the important effects of diurnal weather changes due to cloud and atmospheric scattering were not considered. In this study, we present a SIF upscaling method using photosynthetically active radiation (PAR) as a driving variable. First, a conversion factor (i.e., the ratio of the instantaneous PAR (PARinst) to daily PAR (PARdaily)) was used to upscale in-situ SIF measurements from the instantaneous to daily scale. Then, the performance of the SIF upscaling method was evaluated under changing weather conditions and different latitudes using continuous tower-based measurements at two sites. The results prove that our PAR-based method can reduce not only latitude-dependent but also the weather-dependent variations in the SIF-GPP model. Specifically, the PAR-based method gave a more accurate prediction of diurnal and daily SIF (SIFdaily) than the cos ( SZA ) -based method, with decreased relative root mean square error (RRMSE) values from 42.2% to 25.6% at half-hour intervals and from 25.4% to 13.3% at daily intervals. Moreover, the PAR-based upscaled SIFdaily had a stronger correlation with the daily absorbed PAR (APAR) than both the SIFinst and cos ( SZA ) -based upscaled SIFdaily, especially for cloudy days with a coefficient of determination (R2) that increased from approximately 0.5 to 0.8. Finally, the PAR-based SIFdaily was linked to GPPdaily and compared to the SIFinst or cos ( SZA ) -based SIFdaily. The results indicate that the SIF-GPP correlation can obviously be improved, with an increased R2 from approximately 0.65 to 0.75. Our study confirms the importance of upscaling SIF from the instantaneous to daily scale when linking SIF with GPP and emphasizes the need to take diurnal weather changes into account for SIF temporal upscaling

Atmospheric Correction for Tower-Based Solar-Induced Chlorophyll Fluorescence Observations at O₂-A Band

Solar-induced chlorophyll fluorescence (SIF) has been proven to be an efficient indicator of vegetation photosynthesis. To investigate the relationship between SIF and Gross Primary Productivity (GPP), tower-based continuous spectral observations coordinated with eddy covariance (EC) measurements are needed. As the strong absorption effect at the O2-A absorption bands has an obvious influence on SIF retrieval based on the Fraunhofer Line Discrimination (FLD) principle, atmospheric correction is required even for tower-based SIF observations made with a sensor tens of meters above the canopy. In this study, an operational and simple solution for atmospheric correction of tower-based SIF observations at the O2-A band is proposed. The aerosol optical depth (AOD) and radiative transfer path length (RTPL) are found to be the dominant factors influencing the upward and downward transmittances at the oxygen absorption band. Look-up tables (LUTs) are established to estimate the atmosphere transmittance using AOD and RTPL based on the MODerate resolution atmospheric TRANsmission 5 (MODTRAN 5) model simulations, and the AOD is estimated using the ratio of the downwelling irradiance at 790 nm to that at 660 nm (E790/E660). The influences of the temperature and pressure on the atmospheric transmittance are also compensated for using a corrector factor of RTPL based on an empirical equation. A series of field measurements were carried out to evaluate the performance of the atmospheric correction method for tower-based SIF observations. The difference between the SIF retrieved from tower-based and from ground-based observations decreased obviously after the atmospheric correction. The results indicate that the atmospheric correction method based on a LUT is efficient and also necessary for more accurate tower-based SIF retrieval, especially at the O2-A band

Modeling the Footprint and Equivalent Radiance Transfer Path Length for Tower-Based Hemispherical Observations of Chlorophyll Fluorescence

The measurement of solar-induced chlorophyll fluorescence (SIF) is a new tool for estimating gross primary production (GPP). Continuous tower-based spectral observations together with flux measurements are an efficient way of linking the SIF to the GPP. Compared to conical observations, hemispherical observations made with cosine-corrected foreoptic have a much larger field of view and can better match the footprint of the tower-based flux measurements. However, estimating the equivalent radiation transfer path length (ERTPL) for hemispherical observations is more complex than for conical observations and this is a key problem that needs to be addressed before accurate retrieval of SIF can be made. In this paper, we first modeled the footprint of hemispherical spectral measurements and found that, under convective conditions with light winds, 90% of the total radiation came from an FOV of width 72°, which in turn covered 75.68% of the source area of the flux measurements. In contrast, conical spectral observations covered only 1.93% of the flux footprint. Secondly, using theoretical considerations, we modeled the ERTPL of the hemispherical spectral observations made with cosine-corrected foreoptic and found that the ERTPL was approximately equal to twice the sensor height above the canopy. Finally, the modeled ERTPL was evaluated using a simulated dataset. The ERTPL calculated using the simulated data was about 1.89 times the sensor’s height above the target surface, which was quite close to the results for the modeled ERTPL. Furthermore, the SIF retrieved from atmospherically corrected spectra using the modeled ERTPL fitted well with the reference values, giving a relative root mean square error of 18.22%. These results show that the modeled ERTPL was reasonable and that this method is applicable to tower-based hemispherical observations of SIF

Evaluating the Performance of the SCOPE Model in Simulating Canopy Solar-Induced Chlorophyll Fluorescence

The SCOPE (soil canopy observation of photochemistry and energy fluxes) model has been widely used to interpret solar-induced chlorophyll fluorescence (SIF) and investigate the SIF-photosynthesis links at different temporal and spatial scales in recent years. In the SCOPE model, the fluorescence quantum efficiency in dark-adapted conditions (FQE) for Photosystem II (fqe2) and Photosystem I (fqe1) were two key parameters of SIF emission, which have always been parameterized as fixed values derived from laboratory measurements. To date, only a few studies have focused on evaluating the SCOPE model for SIF interpretation, and the variation of FQE values in the field remains controversial. In this study, the accuracy of the SCOPE model to simulate the canopy SIF was investigated using diurnal experiments on winter wheat. First, ten diurnal experiments were conducted on winter wheat, and the canopy SIF emissions and the SCOPE model’s input parameters were directly measured or indirectly retrieved from the spectral radiances, gross primary productivity (GPP) data, and meteorological records. Second, the SCOPE-simulated SIF emissions with fixed FQE values were evaluated using the observed canopy SIF data. The results show that the SCOPE model can reliably interpret the diurnal cycles of SIF variation and provide acceptable results of SIF simulations at the O2-B (SIFB) and O2-A (SIFA) bands with RRMSEs of 24.35% and 23.67%, respectively. However, the SCOPE-simulated SIFB and SIFA still contained large systematical deviations at some growth stages of wheat, and the seasonal cycles of the ratio between SIFB and SIFA (SIFA/SIFB) cannot be credibly reproduced. Finally, the SCOPE-simulated SIF emissions with variable FQE values were evaluated using the observed canopy SIF data. The simulating accuracy of SIFB and SIFA can be improved greatly using variable FQE values, and the SCOPE simulations track well with the seasonal SIFA/SIFB values with an RRMSE of 20.63%. The results indicated a clear seasonal pattern of FQE values for unbiased SIF simulation: from the erecting to the flowering stage of wheat, the ratio of fqe1 to fqe2 (fqe1/fqe2) gradually increased from 0.05–0.1 to 0.3–0.5, while the fqe2 value decreased from 0.013 to 0.007. Our quantitative results of the model assessment and the FQE adjustment support the use of the SCOPE model as a powerful tool for interpreting the SIF emissions and can serve as a significant reference for future applications of the SCOPE model

Response of Canopy Solar-Induced Chlorophyll Fluorescence to the Absorbed Photosynthetically Active Radiation Absorbed by Chlorophyll

Solar-induced chlorophyll fluorescence (SIF), which can be used as a novel proxy for estimating gross primary production (GPP), can be effectively retrieved using ground-based, airborne and satellite measurements. Absorbed photosynthetically active radiation (APAR) is the key bridge linking SIF and GPP. Remotely sensed SIF at the canopy level ( S I F c a n o p y ) is only a part of the total SIF emission at the photosystem level. An SIF-based model for GPP estimation would be strongly influenced by the fraction of SIF photons escaping from the canopy ( f e s c ). Understanding the response of S I F c a n o p y to the absorbed photosynthetically active radiation absorbed by chlorophyll ( A P A R c h l ) is a key step in estimating GPP but, as yet, this has not been well explored. In this study, we aim to investigate the relationship between remotely sensed S I F c a n o p y and A P A R c h l based on simulations made by the Soil Canopy Observation Photosynthesis Energy fluxes (SCOPE) model and field measurements. First, the ratio of the fraction of the absorbed photosynthetically active radiation absorbed by chlorophyll ( fPAR c h l ) to the fraction of absorbed photosynthetically active radiation absorbed by green leaves ( fPAR g r e e n ) is investigated using a dataset simulated by the SCOPE model. The results give a mean value of 0.722 for Cab at 5 μg cm−2, 0.761 for Cab at 10 μg cm−2 and 0.795 for other Cab content (ranging from 0.71 to 0.81). The response of S I F c a n o p y to A P A R c h l is then explored using simulations corresponding to different biochemical and biophysical conditions and it is found that S I F c a n o p y is well correlated with A P A R c h l . At the O2-A band, for a given plant type, the relationship between S I F c a n o p y and A P A R c h l can be approximately expressed by a linear statistical model even for different values of the leaf area index (LAI) and chlorophyll content, whereas the relationship varies with the LAI and chlorophyll content at the O2-B band. Finally, the response of S I F c a n o p y to A P A R c h l for different leaf angle distribution (LAD) functions is investigated using field observations and simulations; the results show that f e s c is larger for a planophile canopy structure. The values of the ratio of S I F c a n o p y to A P A R c h l are 0.0092 ± 0.0020 , 0.0076 ± 0.0036 and 0.0052 ± 0.0004 μm−1 sr−1 for planophile vegetables/crops, planophile grass and spherical winter wheat, respectively, at the O2-A band. At the O2-B band, the ratios are 0.0063 ± 0.0014 , 0.0049 ± 0.0030 and 0.0033 ± 0.0004 μm−1 sr−1, respectively. The values of this ratio derived from observations agree with simulations, giving values of 0.0055 ± 0.0002 and 0.0068 ± 0.0001 μm−1 sr−1 at the O2-A band and 0.0032 ± 0.0002 and 0.0047 ± 0.0001 μm−1 sr−1 at the O2-B band for spherical and planophile canopies, respectively. Therefore, both the simulations and observations confirm that the relationship between S I F c a n o p y and APAR c h l is species-specific and affected by biochemical components and canopy structure, especially at the O2-B band. It is also very important to correct for reabsorption and scattering of the SIF radiative transfer from the photosystem to the canopy level before the remotely sensed S I F c a n o p y is linked to the GPP

Upscaling Solar-Induced Chlorophyll Fluorescence from an Instantaneous to Daily Scale Gives an Improved Estimation of the Gross Primary Productivity

Solar-induced chlorophyll fluorescence (SIF) is closely linked to the photosynthesis of plants and has the potential to estimate gross primary production (GPP) at different temporal and spatial scales. However, remotely sensed SIF at a ground or space level is usually instantaneous, which cannot represent the daily total SIF. The temporal mismatch between instantaneous SIF (SIFinst) and daily GPP (GPPdaily) impacts their correlation across space and time. Previous studies have upscaled SIFinst to the daily scale based on the diurnal cycle in the cosine of the solar zenith angle ( cos ( SZA ) ) to correct the effects of latitude and length of the day on the variations in the SIF-GPP correlation. However, the important effects of diurnal weather changes due to cloud and atmospheric scattering were not considered. In this study, we present a SIF upscaling method using photosynthetically active radiation (PAR) as a driving variable. First, a conversion factor (i.e., the ratio of the instantaneous PAR (PARinst) to daily PAR (PARdaily)) was used to upscale in-situ SIF measurements from the instantaneous to daily scale. Then, the performance of the SIF upscaling method was evaluated under changing weather conditions and different latitudes using continuous tower-based measurements at two sites. The results prove that our PAR-based method can reduce not only latitude-dependent but also the weather-dependent variations in the SIF-GPP model. Specifically, the PAR-based method gave a more accurate prediction of diurnal and daily SIF (SIFdaily) than the cos ( SZA ) -based method, with decreased relative root mean square error (RRMSE) values from 42.2% to 25.6% at half-hour intervals and from 25.4% to 13.3% at daily intervals. Moreover, the PAR-based upscaled SIFdaily had a stronger correlation with the daily absorbed PAR (APAR) than both the SIFinst and cos ( SZA ) -based upscaled SIFdaily, especially for cloudy days with a coefficient of determination (R2) that increased from approximately 0.5 to 0.8. Finally, the PAR-based SIFdaily was linked to GPPdaily and compared to the SIFinst or cos ( SZA ) -based SIFdaily. The results indicate that the SIF-GPP correlation can obviously be improved, with an increased R2 from approximately 0.65 to 0.75. Our study confirms the importance of upscaling SIF from the instantaneous to daily scale when linking SIF with GPP and emphasizes the need to take diurnal weather changes into account for SIF temporal upscaling

A Reconstructed Global Daily Seamless SIF Product at 0.05 Degree Resolution Based on TROPOMI, MODIS and ERA5 Data

Satellite-derived solar-induced chlorophyll fluorescence (SIF) has been proven to be a valuable tool for monitoring vegetation’s photosynthetic activity at regional or global scales. However, the coarse spatiotemporal resolution or discrete space coverage of most satellite SIF datasets hinders their full potential for studying carbon cycle and ecological processes at finer scales. Although the recent TROPOspheric Monitoring Instrument (TROPOMI) partially addresses this issue, the SIF still has drawbacks in spatial insufficiency and spatiotemporal discontinuities when gridded at high spatiotemporal resolutions (e.g., 0.05°, 1-day or 2-day) due to its nonuniform sampling sizes, swath gaps, and clouds contaminations. Here, we generated a new global SIF product with Seamless spatiotemporal coverage at Daily and 0.05° resolutions (SDSIF) during 2018–2020, using the random forest (RF) approach together with TROPOMI SIF, MODIS reflectance and meteorological datasets. We investigated how the model accuracy was affected by selection of explanatory variables and model constraints. Eventually, models were trained and applied for specific continents and months given the similar response of SIF to environmental variables within closer space and time. This strategy achieved better accuracy (R2 = 0.928, RMSE = 0.0597 mW/m2/nm/sr) than one universal model (R2 = 0.913, RMSE = 0.0653 mW/m2/nm/sr) for testing samples. The SDSIF product can well preserve the temporal and spatial characteristics in original TROPOMI SIF with high temporal correlations (mean R2 around 0.750) and low spatial residuals (less than ±0.081 mW/m2/nm/sr) between them two at most regions (80% of global pixels). Compared with the original SIF at five flux sites, SDSIF filled the temporal gaps and was better consistent with tower-based SIF at the daily scale (the mean R2 increased from 0.467 to 0.744. Consequently, it provided more reliable 4-day SIF averages than the original ones from sparse daily observations (e.g., the R2 at Daman site was raised from 0.614 to 0.837), which resulted in a better correlation with 4-day tower-based GPP. Additionally, the global coverage ratio and local spatial details had also been improved by the reconstructed seamless SIF. Our product has advantages in spatiotemporal continuities and details over the original TROPOMI SIF, which will benefit the application of satellite SIF for understanding carbon cycle and ecological processes at finer spatial and temporal scales

A Reconstructed Global Daily Seamless SIF Product at 0.05 Degree Resolution Based on TROPOMI, MODIS and ERA5 Data

Satellite-derived solar-induced chlorophyll fluorescence (SIF) has been proven to be a valuable tool for monitoring vegetation’s photosynthetic activity at regional or global scales. However, the coarse spatiotemporal resolution or discrete space coverage of most satellite SIF datasets hinders their full potential for studying carbon cycle and ecological processes at finer scales. Although the recent TROPOspheric Monitoring Instrument (TROPOMI) partially addresses this issue, the SIF still has drawbacks in spatial insufficiency and spatiotemporal discontinuities when gridded at high spatiotemporal resolutions (e.g., 0.05°, 1-day or 2-day) due to its nonuniform sampling sizes, swath gaps, and clouds contaminations. Here, we generated a new global SIF product with Seamless spatiotemporal coverage at Daily and 0.05° resolutions (SDSIF) during 2018–2020, using the random forest (RF) approach together with TROPOMI SIF, MODIS reflectance and meteorological datasets. We investigated how the model accuracy was affected by selection of explanatory variables and model constraints. Eventually, models were trained and applied for specific continents and months given the similar response of SIF to environmental variables within closer space and time. This strategy achieved better accuracy (R2 = 0.928, RMSE = 0.0597 mW/m2/nm/sr) than one universal model (R2 = 0.913, RMSE = 0.0653 mW/m2/nm/sr) for testing samples. The SDSIF product can well preserve the temporal and spatial characteristics in original TROPOMI SIF with high temporal correlations (mean R2 around 0.750) and low spatial residuals (less than ±0.081 mW/m2/nm/sr) between them two at most regions (80% of global pixels). Compared with the original SIF at five flux sites, SDSIF filled the temporal gaps and was better consistent with tower-based SIF at the daily scale (the mean R2 increased from 0.467 to 0.744. Consequently, it provided more reliable 4-day SIF averages than the original ones from sparse daily observations (e.g., the R2 at Daman site was raised from 0.614 to 0.837), which resulted in a better correlation with 4-day tower-based GPP. Additionally, the global coverage ratio and local spatial details had also been improved by the reconstructed seamless SIF. Our product has advantages in spatiotemporal continuities and details over the original TROPOMI SIF, which will benefit the application of satellite SIF for understanding carbon cycle and ecological processes at finer spatial and temporal scales